Low- Temperature Low-Field Phases of the Pyrochlore Quantum Magnet Tb 2 Ti 2 07 
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By means of ac magnetic-susceptibility measurements, we find evidence for a new magnetic phase 
of Tb2Ti2C>7 below about 140 mK in zero magnetic field. In magnetic fields parallel to [111], this 
phase is characterized by frequency- and amplitude-dependent susceptibility and extremely slow spin 
dynamics. In the zero-temperature limit, it extends to about 67 mT (the internal field Hi nt — 52 mT), 
at which it makes transition to another phase. The field dependence of the susceptibility of this 
second phase, which extends to about 0.60 T (Hi n ± — 0.54 T) in the zero-temperature limit, indicates 
the presence of a weak magnetization plateau below about 50 mK, as has been predicted by a single- 
tetrahedron four-spin model, suggesting that the second phase is a quantum kagome ice. 

PACS numbers: 75.30.Kz, 75.40.Cx, 75.40.Gb, 75.50.Lk 



In rare-earth-titanatc pyrochlorcs, R,2Ti207, trivalent 
rare-earth ions R 3+ with eightfold oxygen coordina- 
tion form a three-dimensional lattice of corner-sharing 
tctrahcdra. Alternatively, these magnetic oxides can 
be viewed as kagome layers of rare-earth spins coupled 
via interspacing triangular lattices of rare-earth spins, 
stacked along the [111] direction. In either view, the 
geometry leads to frustrated nearest-neighbor exchange 
interactions, whose interplay with an anisotropy dictated 
by a local <111> direction, dipolc-dipole interaction, and 
quantum fluctuations in some cases, result in various ex- 
otic magnetic ground states jl[. Dy 2 Ti2 0_7 and H02T12O7 
are spin ices of classical Ising spins [2 5]. By contrast, 
Gd2Ti2C>7 and Er2Ti2C>7 undergo magnetic transitions at 
temperatures of the order of 1 K, but their ground states 
are unconventional, with fluctuations persisting down to 
20 mK [1,0. The most enigmatic are Tb 2 Ti 2 7 @-[l7| 
and Yb 2 Ti 2 7 [H|-[2l|, which arc under intense investi- 
gation. 

In Tb2Ti207, no long-range order has been found by 
muon-spin relaxation (/zSR) down to 70 mK and by neu- 
tron scattering to 50 mK jg, [l(| , two orders of magni- 
tude lower than the absolute value of the crystal-field- 
subtracted Curie- Weiss temperature, about — 13 K or 
—7. OK [22|, l23j, raising the possibility that this magnet 
is a long-sought three-dimensional quantum spin liquid. 
Indeed, inelastic neutron scattering suggests a crossover 
from a thermall y d isordered paramagnet to a spin liquid 
at about 0.4 K [13] . But a sharp peak in specific heat, 
suggestive of long-range ordering, has been observed at 



0.37 K by a semi-adiabatic method [Uj, although no peak 
has been detected by a relaxation method in a different 
sample [12 ]. On the other hand, a muon-spin-rotation 
(/iSR) frequency shift in low magnetic fields, 20 mT and 
60 mT, applied along the [110] direction suggests a tran- 
sition at a lower temperature of about 150 mK, a transi- 
tion which has not been identified [13]. The static mag- 



netic susceptibility of a zero-field-cooled polycrystallinc 
sample in a 1 mT field shows history dependence sug- 
gestive of spin-glass behavior below about 100 mK, with 
an anomaly at about 70 mK interpreted as the spin-glass 
transition from the high-temperature phase 

Numerical diagonalization of a single-tetrahedron four- 
spin model predicts [24| [25| that the zero-field ground 
state of Tb2Ti2C>7 is a quantum spin liquid, dubbed 
a quantum spin ice p3 . [26|, which — with low increas- 
ing field along the [111] direction — gradually turns into 
a partially polarized state akin to the kagome-ice state 
271 |28| of classical spin-ice magnets. At higher fields, 
H > 82 mT, it evolves into a "three-in one-out" state, 
with three spins pointing into and one pointing out of 
every tetrahedron. The partially polarized state, the 
quantum kagome ice, will manifest itself as a magneti- 
zation plateau spanning from 18 mT to 82 mT, a feature 
detectable at and below 20 mK. 

Motivated by this prediction, we have made system- 
atic measurements of the ac magnetic susceptibility of 
Tb 2 Ti 2 07 in zero magnetic field and fields up to 1.5 T, 
applied along [111]. We find evidence for a new phase 
bounded in zero field at 140 mK and, in the limit of zero 
temperature, at about 67 mT. In this phase, Tb2Ti207 
exhibits frequency- and amplitude-dependent suscepti- 
bility and extremely slow spin dynamics, which is ob- 
served as the temperature is swept or the ac frequency is 
changed stepwise. At 16 mK, the susceptibility indicates 
the presence of a weak magnetization plateau, adjacent 
to the low- field phase and extending to 0.59 T — in qual- 
itative, but not quantitative, agreement with the predic- 
tion. 

The method of sample fabrication is similar to that 
of Ref. [29j]. A single crystal of Tb2Ti2C>7 was grown 
by the floating-zone technique at a rate of 2.5 mm/hour 
in a 0.3 MPa oxygen atmosphere. The sample was ap- 
proximately a square cuboid, 4 mm x 2.4 mm x 2.4 mm, 
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FIG. 1. (color online), (a) x' as a function of temperature 
at various magnetic fields up to 49 mT, measured during up- 
ward temperature sweeps at 0.6mK/min. The frequency and 
amplitude of the ac field were 0.21 Hz and 0.94 mT. (b) \ as 
a function of magnetic field up to 0.25 T at various tempera- 
tures. At and above 85 mK, the frequency and amplitude of 
the ac field and the field-sweep rate were 87.6 Hz, 0.94 mT, 
and 4.2mT/min. For lower temperatures, they were reduced 
to 9.6 Hz, 0.20 mT, and 0.625 mT/min. (With the ac-field am- 
plitude of 0.94 mT, the 100 mK and 85 mK data were taken 
in a slightly nonlinear regime, at least at H = 0. See the Sup- 
plemental Material [H|].) (c) The data at and below 70 mK 
over a wider field range. Lines are field-sweep data. Circles 
are equilibrium data at 40 mK. (d) The 16 mK data up to 
1.5 T, arrows indicating anomalies discussed in the text. See 
the Supplemental Material [33] for x" ■ 



cut from the crystal with the long edges parallel to [111]. 
The setup for ac susceptibility measurements has been 
described elsewhere [3CH32| . The ac field was applied 
along the [111] direction, as was the dc field. The ac 
frequencies ranged from 0.21 Hz to 87.6 Hz, and the ac- 
field amplitude, H ac , was 0.20 mT or 0.94 mT. In addi- 
tion, frequencies up to 389.6 Hz were used to measure 
spin relaxation, and amplitudes ranging from 0.048 mT 
to 1.2 mT were used in a nonlinearity study described in 
the Supplemental Material 



Immersing the sample 
in liquid 3 He, which in turn was cooled by a dilution re- 
frigerator, allowed us to make reliable measurements to 
at least 16 mK, lower than in any previous experiment on 
Tb2Ti 2 07- 

Shown in Fig. Q] are the real part, x'i °f the suscepti- 
bility as functions of temperature and magnetic field, for 
temperatures up to 430 mK and fields up to 1.5 T [34|. 
Because of the extremely slow relaxation of the spins, 
a subject discussed later in detail, the temperature was 
raised very slowly upward at a rate of 0.6mK/min in 
temperature sweeps. For field sweeps, the sweep rate at 
temperatures down to 86 mK was 4.2 mT/min, which was 



reduced to 0.625 mT/min at lower temperatures in order 
to ensure a nearly equilibrium condition (see the Supple- 
mental Material [HI). As a result of the extremely slow 
relaxation, most of our data, even those at 0.21 Hz, were 
taken in the adiabatic limit, not the isothermal limit. 

At zero field, a peak appears in x' a t 154 mK. When 
a small magnetic field, up to 49 mT, is applied along 
[111], the peak becomes smaller and moves to lower tem- 
peratures, as shown in Fig. [lja). In the field sweep 
at 16 mK, a round peak occurs at 63mT, moving to 
lower fields as the temperature is raised, as shown in 
Fig- IHb)- The peak vanishes completely between 85 mK 
and 100 mK. Most significantly, a second peak appears 
at 0.59 T (the internal field Hint — 0.53 T), as shown in 
Fig- HJc) [35|- This smaller, broader peak also moves 
to lower fields as the temperature is raised, becoming a 
very weak shoulder at 50 mK and disappearing before the 
temperature reaches 70 mK. This peak is distinct from 
the weak anomaly observed by Legl et al. jTH in the 
derivative of magnetization, dAI/dH, at a higher field of 
1.9 T (Hi n t - 1-3 T) at 70 mK. It is very likely that that 
shoulder-like feature corresponds to the similar one in 
our data at about 1.2 T (H mt ~ LOT), shown in Fig.QId) 
11. 

The positions of the x' peaks are summarized in Fig. [21 
which suggests the existence of two new phases — phase 
I, bounded at about 140 mK in zero field and at about 
70 mT in the limit of zero temperature, and phase II, 
bounded at about 0.60 T in the same limit. As shown 
in the inset, the upper boundary of phase II exists up 
to about 40 mK but, unlike that of phase I, it gradually 
disappears at higher temperatures. This absence of a 
clear phase boundary to the high-temperature disordered 
phase is reminiscent of a classical kagome ice [28| . 

The field dependence of x' a t 16 mK, with two peaks, 
indicates the presence of a weak magnetization plateau in 
phase II, between 63 mT and 0.59 T. This finding lends 
support to the quantum kagome-ice state predicted by 
the single-tetrahedron four-spin model and, by implica- 
tion, to the underlying proposal [24[ that the low-field 
ground state of Tb2Ti2C>7 is a quantum spin ice. The 
temperature at which the second peak becomes evident, 
40 mK, is consistent with this prediction, strongly sug- 
gesting that the energy scale of the model is correct. The 
two critical fields, however, are higher than predicted, by 
a factor of about three and seven, respectively [35|, sug- 
gesting that the field dependence of the energy level of the 
three-in one-out state may be strongly affected by many- 
tetrahedron effects, which are absent from the model. 

Figures 03a) and E)Jb) show the real and imaginary 
parts, x' an d x", °f the susceptibility near the x' peak 
at zero field and at equilibrium, for three frequencies. At 
0.21 Hz, a round peak appears in x' at about 140 mK. 
With increasing frequency, the x' peak becomes wider 
and moves to slightly higher temperatures. A peak also 
appears in x" about 30 mK below the x' peak. Previous 
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FIG. 2. (color online). H-T phase diagram determined from 
the x peak positions. Triangles are from the temperature 
sweeps, with peak positions in downward and upward sweeps 
averaged. Circles are from field sweeps. Red and magenta 
circles: from the data shown in Figs, [lib) and(c); blue circle: 
from 0.21 Hz data, for which the ac amplitude and the field- 
sweep rate were 0.94 mT and 0.51mT/min. Open symbols 
are the same data plotted against Hint- Lines are guides to 
the eye. Inset: contour plot of \' in the H-T plane in and 
near phase II, constructed from field and temperature sweeps 
at 9.6 Hz. 



susceptibility measurements (loj have observed a similar 
frequency dependence, but with a broad peak at 16 Hz 
appearing at about 250 mK, clearly higher than 160 mK 
at 19.6 Hz for our sample. Measurements by other groups 
[u, 14 1 have also found a similar frequency dependence, 
with peaks in %' and x" occurring at temperatures closer 
to ours, for an ac field applied along [001] or [110]. 

The shift of a x' peak to a higher temperature caused 
by increasing the frequency, observed in these and our 
measurements, is reminiscent of a spin-glass transition 
[36|, as is the appearance of a sharper peak in x" ■ More- 
over, we find that the x'-peak shift per decade frequency 
is ATp/(T p log/) =0.075, where T p is the peak tempera- 
ture and / the frequency, consistent with 0.06-0.08 from 
more extensive measurements [FjJ and similar to that 
of a spin-glass transition in an insulator (3F 



However, 

the downward shift of a x' peak with increasing field 
(Fig. H^a)) is not found in a canonical spin glass, where 
application of a field instead moves the peak slightly up- 
ward in temperature [36|. 

The peak position in field sweeps also depends on the 
frequency, moving to lower fields with increasing fre- 
quency, as shown in Fig. [21c). Therefore, the bound- 
ary of phase I drawn in Fig. [5] is frequency-dependent. 
Moreover, x' at zero field increases linearly with increas- 
ing frequency and exhibits slight hysteresis, as shown in 
Fig. E]Jd). This frequency dependence is much weaker 
at and above 140 mK, outside phase I, and is absent at 
0.36 T, well into phase II. 

Although its mechanism is not understood, this pecu- 
liar frequency dependence of x' allows us to accurately 
measure the relaxation time of spins by stepping the 




FIG. 3. (color online). Equilibrium \' (a) and \" (b) near the 
x' peak in zero field, (c) x' a t 40 mK as a function of magnetic 
field at various frequencies. The field was swept upward at 
4.2mT/min. For all panels, the ac-field amplitude H ac was 
0.94 mT. Thus, in panels (a) and (b), the 87.6 Hz data near 
100 mK were taken in a slightly nonlinear regime. Similarly, 
the 87.6 Hz data in panel (c) were from a slightly nonlinear 
regime, at least at H = (see the Supplemental Material [33|). 
(d) x at H = and 0.36 T as a function of frequency at various 
temperatures. H ac was 0.20 mT. For the 40 mK data at H — 0, 
arrows indicate the directions of frequency steps, showing a 
weak hysteresis. Hysteresis is absent from all the other data. 
See the Supplemental Material [33] for the imaginary part of 
the data shown in panel (c). 



frequency, which — unlike temperature — can be changed 
nearly instantaneously. Figure HJa) shows the total mag- 
nitude, |x|, of the susceptibility at 40 mK at zero field, as 
the frequency was changed stepwise between 9.6 Hz and 
389.6 Hz, first upward then downward. instead of x' 
alone, is fitted to a stretched exponential cxp[— (t/r) a ] 
to extract the relaxation time r, because both x' an d x" 
evolve with time after a frequency step. The stretching 
exponent a ranges from 0.94 to 1.10, indicating that the 
relaxation is nearly exponential. The relaxation time is 
a strong function of frequency at this temperature, as 
shown in Fig. 2th). It also depends strongly on tempera- 
ture, as Fig.[4tc) shows, with an abrupt rise within a nar- 
row temperature interval between 120 mK and 130 mK, 
where x" rises sharply (see Fig. EJb)). The relaxation 
time measured at 0.21 Hz up to 120 mK can be expressed 
by an empirical formula for spin glasses |37|, 1/t ~ 
exp[— a(T p /T)] with a = 0.41 and T p = 140mK, as shown 
in Fig. 0|d). Note that even at 140 mK, the relaxation 
rate, 1/t, is as small as 7.1 x 10 -3 s _1 , two orders of mag- 
nitude smaller than even our lowest measuring frequency, 
indicating that much of the susceptibility we have mea- 
sured is in the adiabatic limit, not the isothermal limit. 

A similar temperature dependence has been found in 
a spin glass (38l . [39) . Unlike in our case, however, the 
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FIG. 4. (color online), (a) \\\ as a function of time at 
40 mK in zero field. Before t = 0, the ac frequency was held 
at 2.6 Hz for 12 hours until equilibrium was reached. There- 
after, the frequency was changed stepwise, each time after 
new equilibrium was reached. For clarity, relaxation curves 
at 189.6 Hz, 269.6 Hz, and 389.6 Hz— between the two 87.6 Hz 
curves — are not shown, (b) Relaxation time r as a function 
of frequency, extracted from the data. Arrows indicate the 
frequency-step sequence, (c) r vs temperature. At each tem- 
perature, the frequency was first held at 389.6 Hz until equi- 
librium was reached, then stepped down to 0.21 Hz, the fre- 
quency at which the measurements were made, (d) 1/r vs 
T p /T, where T p = 140 mK. For all the panels, H ac = 0.20 mT. 



relaxation is logarithmic in time, as predicted by the 
droplet scaling theory [40j . not nearly exponential. To 
our knowledge, relaxation rates much smaller than mea- 
suring frequencies have not been observed in a spin glass. 

In summary, we have constructed, from ac magnetic 
susceptibility, a low-temperature magnetic phase dia- 
gram of Tb2Ti 2 07 in low magnetic fields applied along 
[111], uncovering two phases. The field dependence of the 
susceptibility at 16 mK and 40 mK indicates the presence 
of a weak magnetization plateau in support of theory 
24, 25|, suggesting that phase II is a quantum kagome 
ice. In view of the underlying proposal of the theory, 
phase I is likely to be a quantum spin ice. The fre- 
quency dependence of the susceptibility-peak position at 
zero field is similar to that of a spin glass, as is the tem- 
perature dependence of the relaxation time. But the field 
dependence of the peak position and the time dependence 
of the relaxation, as well as the frequency and amplitude 
dependence of x' at 40 mK, suggest that phase I is not a 
spin glass. 

In addition to our work, three very recent studies have 
investigated Tb2Ti 2 C>7 in magnetic fields applied along 
[111], also motivated by the prediction on the low-field 
ground state by the theory. Lhotel et al. 14j measured 
the magnetization up to 8T and ac susceptibilities at 
zero field, and found no evidence for a magnetization 
plateau. However, the lowest temperature of the experi- 



ment was 80 mK, insufficient to detect the weak plateau 
indicated by our susceptibility result. Legl et al. [15| also 
measured the magnetization, to a lower temperature of 
43 mK, yet finding no evidence for a plateau. In par- 
ticular, the derivative dM/dH exhibits no anomaly that 
corresponds to the second \' peak we have observed at 
and below 40 mK. Our results shown in Fig.QJc) suggest, 
however, that 43 mK may not be cold enough to clearly 
detect such an anomaly. Moreover, the field-sweep rate 
employed by Legl et al. was four times faster than the 
rate we have found at 40 mK to completely wipe out the 
peak (see the Supplemental Material [33| )• The low-field 
low-temperature phase diagram proposed by these au- 
thors also differs substantially from ours. However, at 
least part of the discrepancy may come from the much 
faster temperature-sweep rate — 5 mK/min as opposed to 
our 0.6 mK/min — used by Legl et al. Baker et al. [l6| 
measured //SR to 25 mK and ac susceptibility to 68 mK. 
They have observed anomalies in the /iSR relaxation time 
at 15 mT and ~60mT, at and below 50 mK. The sus- 
ceptibility x' a t 68 mK and 100 mK, measured at 50 Hz 
and 500 Hz, is similar to ours but shows a broad peak 
at about 15mT, lower than the positions of the peak we 
have found at 9.6 Hz and 87.6 Hz at similar temperatures 
(Fig. Gib)). The temperature and field ranges of the ex- 
periment were not sufficient to observe the second peak, 
the key signature of the weak magnetization plateau. 
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Effect of slow relaxation. Spin relaxation in Tb2Ti207 
slows down dramatically in phase I as the temperature 
decreases, as shown in Fig. 4(c) of the accompanying Let- 
ter. Therefore, great care must be exercised to ensure 
equilibrium when measurements are made while sweep- 
ing the temperature or the magnetic field. 
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FIG. SI: (color online), (a) Real part, \ , and (b) imaginary 
part, x"t °f the ac susceptibility at 0.21 Hz as a function of 
temperature at zero field, with arrows indicating the direc- 
tions of temperature sweeps. Blue circles are the equilibrium 
data shown in Figs. 3(a) and 3(b). The ac-field amplitude 
was # ac = 0.94mT. 

The effect of the slow relaxation on ac susceptibility 
measured during temperature sweeps was studied at zero 
field with an ac frequency of 0.21 Hz and an amplitude 
of 0.94 mT. Starting from 500 mK, the temperature was 
lowered to 25 mK, where it was held for eight hours, then 
raised back to 500 mK. The temperature regulator was 
set at a sweep rate of 0.6 mK/min, but the actual cooling 
rate gradually decreased below 100 mK, limited by the 



cooling power of the dilution refrigerator; for instance, 
at 70 mK it was 0.2 mK/min. 

Below about 150 mK, the cooling and warming curves 
of both x' and x" 01 the ac susceptibility deviate strongly 
from each other, as shown in Fig. [Si] [1(. But the equi- 
librium data taken near the x' peak shows no hysteresis 
within our resolution, indicating that the hysteresis arises 
from slow spin relaxation. Indeed, at the lowest temper- 
ature, 25 mK, the relaxation time of x' was as long as 
6.8 hours. Because the warming rate was constant, un- 
like the cooling rate, only data taken upon warming arc 
shown in Fig. 1(a) of the accompanying Letter. The data 
shown in Fig. [Si] suggest, however, that the average of 
the x'-peak positions in the cooling and warming curves 
is very close to that of the equilibrium data. In Fig. 2 
of the Letter, such averages are used in part to construct 
the phase diagram. 
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FIG. S2: (color online). x' as a function of magnetic 
field at 40 mK, for two field-sweep rates in comparison with 
equilibrium data. The magnetic field, applied along the 
[ill] direction, was swept upward. The ac frequency was 
9.6 Hz, and H ac was 0.20 mT. The equilibrium data and the 
0.625 mT/min data are shown also in Fig. 1(c). 

Great care must be exercised also when data are taken 
while sweeping the field. The effect of the slow relaxation 
on field-sweep data was studied at 40 mK with an ac fre- 
quency of 9.6 Hz and an amplitude of 0.20 mT. We have 



2 



found that even a sweep rate as low as 3.75mT/min, a 
quarter of the rate used by Legl. et al. Q, completely 
wipes out the 0.52 T peak in x'> whereas a field sweep at 
0.625 mT/min reveals the peak in agreement with equi- 
librium data taken at fixed fields, as shown in Fig. [S2] 
For this reason, all the field-sweep data shown in Fig. 1(c) 
of the accompanying Letter were taken with a sweep rate 
of 0.625 mT/min. 

Temperature sweeps at and above 0.15 T. Figure [S3] 
shows x' an d x" m the field region of phase II and, for 
comparison, at 0.8 T, which is slightly above phase II. 
No anomaly appears in these temperature-sweep data, 
indicating that the boundary of phase II is open, unlike 
the closed boundary of phase I. x" is very small at these 
fields; as a result, it sometimes picks up semi-periodic 
noise, which may come from the air conditioning of the 
laboratory. 
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FIG. S3: (color online), x an d x" as a function of tempera- 
ture at various fields. The frequency and amplitude of the ac 
field were 9.6 Hz and 0.20 mT. The temperature was raised at 
a rate of 0.7mK/min. 

Amplitude dependence of susceptibility in phase I. In 
zero field, and at and blow 100 mK, x' depends on the ac- 
field amplitude H ac . The amplitude dependence has an 
onset, as shown in Fig.lS4l At 87.6 Hz, the onset occurs at 
if ac ~0.9mT, above which x' increases with increasing 
amplitude, whereas the dependence is weaker at 9.6 Hz, 
with the onset moving to _ff ac c±l.lmT. At and above 



150 mK, outside phase I, we do not see an amplitude 
dependence even at 87.6 Hz, at least up to H ac = 1.2 mT. 
The mechanism of this peculiar H ac dependence is not 
understood. 

All the data shown in the accompanying Letter were 
taken in the amplitude-independent, linear regime — 
except for the 100 mK and 85 mK data in Fig. 1(b), the 
87.6 Hz data near 100 mK in Figs. 3(a) and 3(b), and the 
87.6 Hz data in Fig. 3(c). Those data were taken with 
H ac = 0.94 mT, which is in a slightly nonlinear regime at 
least at H = 0. 
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FIG. S4: (color online). Dependence of %' at H — on the 
ac-field amplitude, H ac , at 87.6 Hz (solid circles) and 9.6 Hz 
(open circles). 



Imaginary part of susceptibility. For completeness, we 
present the imaginary part, x"i of the susceptibility data 
whose real part is shown in Figs. 1(a), 1(b), 1(d), and 3(c) 
of the accompanying Letter. As shown in Fig. IS5fa). the 
peak in x" is sharper than the x' peak (see Fig. 1(a)) and 
appears about 28 mK below the latter. With increasing 
field, the x" peak moves to lower temperatures along with 
the x' peak, while at the same time becoming smaller. 
In field sweeps, the x" peak is roughly as broad as the x' 
peak, appearing about 6 mT above the latter. It moves 
to lower fields along with the latter, as the temperature 
is raised, as shown in Fig. lSBT b). x" is negligible at fields 
above about 0.2 T. As shown in Fig. IS61 the effect of 
raising the ac frequency is qualitatively similar to that of 
raising the temperature; with increasing frequency, the 
broad x" peak moves to lower fields along with the x' 
peak. 
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FIG. S6: (color online). Imaginary part, x", of the ac suscep- 
tibility whose real part is shown in Fig. 3(c). 
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All x' an d x" m the Supplemental Material and the accom- 
panying Letter have been corrected for the demagnetizing 
effect. 
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FIG. S5: (color online), (a) Imaginary part, x"> °f the ac 
susceptibility whose real part is shown in Fig. 1(a) of the 
accompanying Letter, (b) x" of the susceptibility whose real 
part is shown in Figs. 1(b) and 1(d). The small hump at a 
field below 0.02 T comes from the background without the 
sample. 



